Dependence of nucleosome mechanical stability on DNA mismatches.

The organization of nucleosomes into chromatin and their accessibility are shaped by local DNA mechanics. Conversely, nucleosome positions shape genetic variations, which may originate from mismatches during replication and chemical modification of DNA. To investigate how DNA mismatches affect the mechanical stability and the exposure of nucleosomal DNA, we used an optical trap combined with single-molecule FRET and a single-molecule FRET cyclization assay. We found that a single base-pair C-C mismatch enhances DNA bendability and nucleosome mechanical stability for the 601-nucleosome positioning sequence. An increase in force required for DNA unwrapping from the histone core is observed for single base-pair C-C mismatches placed at three tested positions: at the inner turn, at the outer turn, or at the junction of the inner and outer turn of the nucleosome. The results support a model where nucleosomal DNA accessibility is reduced by mismatches, potentially explaining the preferred accumulation of single-nucleotide substitutions in the nucleosome core and serving as the source of genetic variation during evolution and cancer progression. Mechanical stability of an intact nucleosome, that is mismatch-free, is also dependent on the species as we find that yeast nucleosomes are mechanically less stable and more symmetrical in the outer turn unwrapping compared to Xenopus nucleosomes.

Activation of coagulation FXI promotes endothelial inflammation and amplifies platelet activation in a nonhuman primate model of hyperlipidemia.

Background: Hyperlipidemia is associated with chronic inflammation and thromboinflammation. This is an underlying cause of several cardiovascular diseases, including atherosclerosis. In diseased blood vessels, rampant thrombin generation results in the initiation of the coagulation cascade, activation of platelets, and endothelial cell dysfunction. Coagulation factor (F) XI represents a promising therapeutic target to reduce thromboinflammation, as it is uniquely positioned at an intersection between inflammation and thrombin generation. Objectives: This study aimed to investigate the role of FXI in promoting platelet and endothelial cell activation in a model of hyperlipidemia. Methods: Nonhuman primates (NHPs) were fed a standard chow diet (lean, n = 6) or a high-fat diet (obese, n = 8) to establish a model of hyperlipidemia. Obese NHPs were intravenously administered a FXI blocking antibody (2 mg/kg) and studied at baseline and at 1, 7, 14, 21, and 28 days after drug administration. Platelet activation and inflammatory markers were measured using fluorescence-activated cell sorting or enzyme-linked immunosorbent assay. Molecular imaging was used to quantify vascular cell adhesion molecule 1 (VCAM-1) expression at the carotid bifurcation. Results: Obese NHPs demonstrated increased sensitivity for platelet P-selectin expression and phosphatidylserine exposure in response to platelet GPVI or PAR agonists compared with lean NHPs. Obese NHPs exhibited elevated levels of C-reactive protein, cathepsin D, and myeloperoxidase compared with lean NHPs. Following pharmacological inhibition of FIX activation by FXIa, platelet priming for activation by GPVI or PAR agonists, C-reactive protein levels, and endothelial VCAM-1 levels were reduced in obese NHPs. Conclusion: FXI activation promotes the proinflammatory phenotype of hyperlipidemia by priming platelet activation and inciting endothelial cell dysfunction.

Selective enrichment of plasma cell-free messenger RNA in cancer-associated extracellular vesicles.

Extracellular vesicles (EVs) have been shown as key mediators of extracellular small RNA transport. However, carriers of cell-free messenger RNA (cf-mRNA) in human biofluids and their association with cancer remain poorly understood. Here, we performed a transcriptomic analysis of size-fractionated plasma from lung cancer, liver cancer, multiple myeloma, and healthy donors. Morphology and size distribution analysis showed the successful separation of large and medium particles from other soluble plasma protein fractions. We developed a strategy to purify and sequence ultra-low amounts of cf-mRNA from particle and protein enriched subpopulations with the implementation of RNA spike-ins to control for technical variability and to normalize for intrinsic drastic differences in cf-mRNA amount carried in each plasma fraction. We found that the majority of cf-mRNA was enriched and protected in EVs with remarkable stability in RNase-rich environments. We observed specific enrichment patterns of cancer-associated cf-mRNA in each particle and protein enriched subpopulation. The EV-enriched differentiating genes were associated with specific biological pathways, such as immune systems, liver function, and toxic substance regulation in lung cancer, liver cancer, and multiple myeloma, respectively. Our results suggest that dissecting the complexity of EV subpopulations illuminates their biological significance and offers a promising liquid biopsy approach.

Plasma cell-free RNA profiling distinguishes cancers from pre-malignant conditions in solid and hematologic malignancies.

Cell-free RNA (cfRNA) in plasma reflects phenotypic alterations of both localized sites of cancer and the systemic host response. Here we report that cfRNA sequencing enables the discovery of messenger RNA (mRNA) biomarkers in plasma with the tissue of origin-specific to cancer types and precancerous conditions in both solid and hematologic malignancies. To explore the diagnostic potential of total cfRNA from blood, we sequenced plasma samples of eight hepatocellular carcinoma (HCC) and ten multiple myeloma (MM) patients, 12 patients of their respective precancerous conditions, and 20 non-cancer (NC) donors. We identified distinct gene sets and built classification models using Random Forest and linear discriminant analysis algorithms that could distinguish cancer patients from premalignant conditions and NC individuals with high accuracy. Plasma cfRNA biomarkers of HCC are liver-specific genes and biomarkers of MM are highly expressed in the bone marrow compared to other tissues and are related to cell cycle processes. The cfRNA level of these biomarkers displayed a gradual transition from noncancerous states through precancerous conditions and cancer. Sequencing data were cross-validated by quantitative reverse transcription PCR and cfRNA biomarkers were validated in an independent sample set (20 HCC, 9 MM, and 10 NC) with AUC greater than 0.86. cfRNA results observed in precancerous conditions require further validation. This work demonstrates a proof of principle for using mRNA transcripts in plasma with a small panel of genes to distinguish between cancers, noncancerous states, and precancerous conditions.

Irreversible alteration of extracellular vesicle and cell-free messenger RNA profiles in human plasma associated with blood processing and storage.

The discovery and utility of clinically relevant circulating biomarkers depend on standardized methods that minimize preanalytical errors. Despite growing interest in studying extracellular vesicles (EVs) and cell-free messenger RNA (cf-mRNA) as potential biomarkers, how blood processing and freeze/thaw impacts the profiles of these analytes in plasma was not thoroughly understood. We utilized flow cytometric analysis to examine the effect of differential centrifugation and a freeze/thaw cycle on EV profiles. Utilizing flow cytometry postacquisition analysis software (FCMpass) to calibrate light scattering and fluorescence, we revealed how differential centrifugation and post-freeze/thaw processing removes and retains EV subpopulations. Additionally, cf-mRNA levels measured by RT-qPCR profiles from a panel of housekeeping, platelet, and tissue-specific genes were preferentially affected by differential centrifugation and post-freeze/thaw processing. Critically, freezing plasma containing residual platelets yielded irreversible ex vivo generation of EV subpopulations and cf-mRNA transcripts, which were not removable by additional processing after freeze/thaw. Our findings suggest the importance of minimizing confounding variation attributed to plasma processing and platelet contamination.

Longitudinal Analysis of Human Pancreatic Adenocarcinoma Development Reveals Transient Gene Expression Signatures.

Previous transcriptome studies of human pancreatic ductal adenocarcinoma (PDAC) compare non-cancerous pancreatic intraepithelial neoplasias (PanIN) with late-stage PDAC obtained from different patients, thus have limited ability to discern network dynamics that contribute to the disease progression. We demonstrated previously that the 10-22 cell line, an induced pluripotent stem cell-like line reprogrammed from late-stage human PDAC cells, recapitulated the progression from PanINs to PDAC upon transplantation into NOD/LtSz-scid/IL2R-gammanull mice. Herein, we investigated the transition from precursor to PDAC using the isogenic model. We analyzed transcriptomes of genetically tagged 10-22 cells progressing from PanINs to PDAC in mice and validated the results using The Cancer Genome Atlas PDAC dataset, human clinical PanIN and PDAC tissues, and a well-established murine PDAC model. We functionally studied candidate proteins using human normal (H6C7) and cancerous (Miapaca2, Aspc1) pancreatic ductal epithelial cell lines. 10-22 cell-derived PDAC displayed the molecular signature of clinical human PDAC. Expression changes of many genes were transient during PDAC progression. Pathways for extracellular vesicle transport and neuronal cell differentiation were derepressed in the progression of PanINs to PDAC. HMG-box transcription factor 1 (HBP1) and BTB domain and CNC homolog 1 (BACH1) were implicated in regulating dynamically expressed genes during PDAC progression, and their expressions inversely correlated with PDAC patients' prognosis. Ectopic expression of HBP1 increased proliferation and migration of normal and cancerous pancreatic cells, indicating that HBP1 may confer the cell dissemination capacity in early PDAC progression. This unique longitudinal analysis provides insights into networks underlying human PDAC progression and pathogenesis. IMPLICATIONS: Manipulation of HBP1, BACH1, and RUN3 networks during PDAC progression can be harnessed to develop new targets for treating PDAC.

Measuring DNA mechanics on the genome scale.

Mechanical deformations of DNA such as bending are ubiquitous and have been implicated in diverse cellular functions1. However, the lack of high-throughput tools to measure the mechanical properties of DNA has limited our understanding of how DNA mechanics influence chromatin transactions across the genome. Here we develop 'loop-seq'-a high-throughput assay to measure the propensity for DNA looping-and determine the intrinsic cyclizabilities of 270,806 50-base-pair DNA fragments that span Saccharomyces cerevisiae chromosome V, other genomic regions, and random sequences. We found sequence-encoded regions of unusually low bendability within nucleosome-depleted regions upstream of transcription start sites (TSSs). Low bendability of linker DNA inhibits nucleosome sliding into the linker by the chromatin remodeller INO80, which explains how INO80 can define nucleosome-depleted regions in the absence of other factors2. Chromosome-wide, nucleosomes were characterized by high DNA bendability near dyads and low bendability near linkers. This contrast increases for deeper gene-body nucleosomes but disappears after random substitution of synonymous codons, which suggests that the evolution of codon choice has been influenced by DNA mechanics around gene-body nucleosomes. Furthermore, we show that local DNA mechanics affect transcription through TSS-proximal nucleosomes. Overall, this genome-scale map of DNA mechanics indicates a 'mechanical code' with broad functional implications.

Diurnal stability of cell-free DNA and cell-free RNA in human plasma samples.

Many emerging technologies are reliant on circulating cell-free DNA (cfDNA) and cell-free RNA (cfRNA) applications in the clinic. However, the impact of diurnal cycles or daily meals on circulating analytes are poorly understood and may be confounding factors when developing diagnostic platforms. To begin addressing this knowledge gap, we obtained plasma from four healthy donors serially sampled five times during 12 h in a single day. For all samples, we measured concentrations of cfDNA and cfRNA using both bulk measurements and gene-specific digital droplet PCR. We found no significant variation attributed to blood draw number for the cfDNA or cfRNA. This indicated that natural diurnal cycles and meal consumption do not appear to significantly affect abundance of total cfDNA, total cfRNA, or our two selected cfRNA transcripts. Conversely, we observed significant variation between individual donors for cfDNA and one of the cfRNA transcripts. The results of this work suggest that it will be important to consider patient-specific baselines when designing reliable circulating cfDNA or cfRNA clinical assays.

Extreme mechanical diversity of human telomeric DNA revealed by fluorescence-force spectroscopy.

G-quadruplexes (GQs) can adopt diverse structures and are functionally implicated in transcription, replication, translation, and maintenance of telomere. Their conformational diversity under physiological levels of mechanical stress, however, is poorly understood. We used single-molecule fluorescence-force spectroscopy that combines fluorescence resonance energy transfer with optical tweezers to measure human telomeric sequences under tension. Abrupt GQ unfolding with K+ in solution occurred at as many as four discrete levels of force. Added to an ultrastable state and a gradually unfolding state, there were six mechanically distinct structures. Extreme mechanical diversity was also observed with Na+, although GQs were mechanically weaker. Our ability to detect small conformational changes at low forces enabled the determination of refolding forces of about 2 pN. Refolding was rapid and stochastically redistributed molecules to mechanically distinct states. A single guanine-to-thymine substitution mutant required much higher ion concentrations to display GQ-like unfolding and refolded via intermediates, contrary to the wild type. Contradicting an earlier proposal, truncation to three hexanucleotide repeats resulted in a single-stranded DNA-like mechanical behavior under all conditions, indicating that at least four repeats are required to form mechanically stable structures.

Multiomics modeling of the immunome, transcriptome, microbiome, proteome and metabolome adaptations during human pregnancy.

Motivation: Multiple biological clocks govern a healthy pregnancy. These biological mechanisms produce immunologic, metabolomic, proteomic, genomic and microbiomic adaptations during the course of pregnancy. Modeling the chronology of these adaptations during full-term pregnancy provides the frameworks for future studies examining deviations implicated in pregnancy-related pathologies including preterm birth and preeclampsia. Results: We performed a multiomics analysis of 51 samples from 17 pregnant women, delivering at term. The datasets included measurements from the immunome, transcriptome, microbiome, proteome and metabolome of samples obtained simultaneously from the same patients. Multivariate predictive modeling using the Elastic Net (EN) algorithm was used to measure the ability of each dataset to predict gestational age. Using stacked generalization, these datasets were combined into a single model. This model not only significantly increased predictive power by combining all datasets, but also revealed novel interactions between different biological modalities. Future work includes expansion of the cohort to preterm-enriched populations and in vivo analysis of immune-modulating interventions based on the mechanisms identified. Availability and implementation: Datasets and scripts for reproduction of results are available through: https://nalab.stanford.edu/multiomics-pregnancy/. Supplementary information: Supplementary data are available at Bioinformatics online.

Noninvasive blood tests for fetal development predict gestational age and preterm delivery.

Noninvasive blood tests that provide information about fetal development and gestational age could potentially improve prenatal care. Ultrasound, the current gold standard, is not always affordable in low-resource settings and does not predict spontaneous preterm birth, a leading cause of infant death. In a pilot study of 31 healthy pregnant women, we found that measurement of nine cell-free RNA (cfRNA) transcripts in maternal blood predicted gestational age with comparable accuracy to ultrasound but at substantially lower cost. In a related study of 38 women (23 full-term and 15 preterm deliveries), all at elevated risk of delivering preterm, we identified seven cfRNA transcripts that accurately classified women who delivered preterm up to 2 months in advance of labor. These tests hold promise for prenatal care in both the developed and developing worlds, although they require validation in larger, blinded clinical trials.

Simultaneously Monitoring Immune Response and Microbial Infections during Pregnancy through Plasma cfRNA Sequencing.

BACKGROUND: Plasma cell-free RNA (cfRNA) encompasses a broad spectrum of RNA species that can be derived from both human cells and microbes. Because cfRNA is fragmented and of low concentration, it has been challenging to profile its transcriptome using standard RNA-seq methods. METHODS: We assessed several recently developed RNA-seq methods on cfRNA samples. We then analyzed the dynamic changes of both the human transcriptome and the microbiome of plasma during pregnancy from 60 women. RESULTS: cfRNA reflects a well-orchestrated immune modulation during pregnancy: an up-regulation of antiinflammatory genes and an increased abundance of antimicrobial genes. We observed that the plasma microbiome remained relatively stable during pregnancy. The bacteria Ureaplasma shows an increased prevalence and increased abundance at postpartum, which is likely to be associated with postpartum infection. We demonstrated that cfRNA-seq can be used to monitor viral infections. We detected a number of human pathogens in our patients, including an undiagnosed patient with a high load of human parvovirus B19 virus (B19V), which is known to be a potential cause of complications in pregnancy. CONCLUSIONS: Plasma cfRNA-seq demonstrates the potential to simultaneously monitor immune response and microbial infections during pregnancy.

Defining Single Molecular Forces Required for Notch Activation Using Nano Yoyo.

Notch signaling, involved in development and tissue homeostasis, is activated at the cell-cell interface through ligand-receptor interactions. Previous studies have implicated mechanical forces in the activation of Notch receptor upon binding to its ligand. Here we aimed to determine the single molecular force required for Notch activation by developing a novel low tension gauge tether (LTGT). LTGT utilizes the low unbinding force between single-stranded DNA (ssDNA) and Escherichia coli ssDNA binding protein (SSB) (∼4 pN dissociation force at 500 nm/s pulling rate). The ssDNA wraps around SSB and, upon application of force, unspools from SSB, much like the unspooling of a yoyo. One end of this nano yoyo is attached to the surface though SSB, while the other end presents a ligand. A Notch receptor, upon binding to its ligand, is believed to undergo force-induced conformational changes required for activating downstream signaling. If the required force for such activation is larger than 4 pN, ssDNA will unspool from SSB, and downstream signaling will not be activated. Using these LTGTs, in combination with the previously reported TGTs that rupture double-stranded DNA at defined forces, we demonstrate that Notch activation requires forces between 4 and 12 pN, assuming an in vivo loading rate of 60 pN/s. Taken together, our study provides a direct link between single-molecular forces and Notch activation.

Effects of cytosine modifications on DNA flexibility and nucleosome mechanical stability.

Cytosine can undergo modifications, forming 5-methylcytosine (5-mC) and its oxidized products 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). Despite their importance as epigenetic markers and as central players in cellular processes, it is not well understood how these modifications influence physical properties of DNA and chromatin. Here we report a comprehensive survey of the effect of cytosine modifications on DNA flexibility. We find that even a single copy of 5-fC increases DNA flexibility markedly. 5-mC reduces and 5-hmC enhances flexibility, and 5-caC does not have a measurable effect. Molecular dynamics simulations show that these modifications promote or dampen structural fluctuations, likely through competing effects of base polarity and steric hindrance, without changing the average structure. The increase in DNA flexibility increases the mechanical stability of the nucleosome and vice versa, suggesting a gene regulation mechanism where cytosine modifications change the accessibility of nucleosomal DNA through their effects on DNA flexibility.

Nucleosomes undergo slow spontaneous gaping.

In eukaryotes, DNA is packaged into a basic unit, the nucleosome which consists of 147 bp of DNA wrapped around a histone octamer composed of two copies each of the histones H2A, H2B, H3 and H4. Nucleosome structures are diverse not only by histone variants, histone modifications, histone composition but also through accommodating different conformational states such as DNA breathing and dimer splitting. Variation in nucleosome structures allows it to perform a variety of cellular functions. Here, we identified a novel spontaneous conformational switching of nucleosomes under physiological conditions using single-molecule FRET. Using FRET probes placed at various positions on the nucleosomal DNA to monitor conformation of the nucleosome over a long period of time (30-60 min) at various ionic conditions, we identified conformational changes we refer to as nucleosome gaping. Gaping transitions are distinct from nucleosome breathing, sliding or tightening. Gaping modes switch along the direction normal to the DNA plane through about 5-10 angstroms and at minutes (1-10 min) time scale. This conformational transition, which has not been observed previously, may be potentially important for enzymatic reactions/transactions on nucleosomal substrate and the formation of multiple compression forms of chromatin fibers.

Asymmetric unwrapping of nucleosomes under tension directed by DNA local flexibility.

Dynamics of the nucleosome and exposure of nucleosomal DNA play key roles in many nuclear processes, but local dynamics of the nucleosome and its modulation by DNA sequence are poorly understood. Using single-molecule assays, we observed that the nucleosome can unwrap asymmetrically and directionally under force. The relative DNA flexibility of the inner quarters of nucleosomal DNA controls the unwrapping direction such that the nucleosome unwraps from the stiffer side. If the DNA flexibility is similar on two sides, it stochastically unwraps from either side. The two ends of the nucleosome are orchestrated such that the opening of one end helps to stabilize the other end, providing a mechanism to amplify even small differences in flexibility to a large asymmetry in nucleosome stability. Our discovery of DNA flexibility as a critical factor for nucleosome dynamics and mechanical stability suggests a novel mechanism of gene regulation by DNA sequence and modifications.

A Coarse-Grained Model of Unstructured Single-Stranded DNA Derived from Atomistic Simulation and Single-Molecule Experiment.

A simple coarse-grained model of single-stranded DNA (ssDNA) was developed, featuring only two sites per nucleotide that represent the centers of mass of the backbone and sugar/base groups. In the model, the interactions between sites are described using tabulated bonded potentials optimized to reproduce the solution structure of DNA observed in atomistic molecular dynamics simulations. Isotropic potentials describe nonbonded interactions, implicitly taking into account the solvent conditions to match the experimentally determined radius of gyration of ssDNA. The model reproduces experimentally measured force-extension dependence of an unstructured DNA strand across 2 orders of magnitude of the applied force. The accuracy of the model was confirmed by measuring the end-to-end distance of a dT14 fragment via FRET while stretching the molecules using optical tweezers. The model offers straightforward generalization to systems containing double-stranded DNA and DNA binding proteins.

Molecular restructuring of water and lipids upon the interaction of DNA with lipid monolayers.

Understanding the molecular mechanism of DNA/lipid interaction is critical in optimizing the use of lipid cofactors in gene therapy. Here, we address this question by employing label-free vibrational sum frequency (VSF) spectroscopy to study the interaction of DNA with lipid monolayers of the cationic lipids DPTAP(1,2-dipalmitoyl-3-trimethylammonium-propane) and diC14-amidine as well as the zwitterionic lipid DPPC (1,2-dipalmitoyl-sn-glycero-3-phosphocholine) in the presence and absence of calcium. Our approach has the advantage both of allowing us to explicitly probe intermolecular interactions and of providing insight into the structure of water and lipids around DNA at the lipid interface. We find, by examination of the OD stretch of interfacial D(2)O, that water structure differs markedly between systems containing DNA adsorbed to cationic and those that contain DNA adsorbed to zwitterionic lipid monolayers (in the presence or absence of Ca(2+)). The spectral response of interfacial water in the cationic system is consistent with a highly structured, undercoordinated, structural 'type' of water. Further, by investigation of CH stretch modes of the diC14-amidine lipid tails, we demonstrate that the adsorption of DNA to this lipid leads to increased ordering of lipid tails.